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Abstract 

We review some of the history and properties of theories for the vari- 
ation of the gravitation and fine structure 'constants'. We highlight some 
general features of the cosmological models that exist in these theories 
with reference to recent quasar data that is consistent with time- variation 
in alpha since a redshift of 3.5. The behaviour of a simple class of varying- 
alpha cosmologies is outlined in the light of all the observational con- 
straints. We discuss the key role played by non-zero vacuum energy and 
curvature in turning off the variation of constants in these theories and the 
issue of comparing extra-galactic and local observational data. We also 
show why black hole thermodynamics does not enable us to distinguish 
between time variations of different constants. 

Introduction 



There are a number of reasons why the possibility of varying constants should 
be taken seriously First, we know that the best candidates for unification 
of the forces of nature in a quantum gravitational environment only seem to 
exist in finite form if there are many more dimensions of space than the three 
that we are familiar with. This means that the true constants of nature are 
defined in higher dimensions and the three-dimensional shadows we observe are 
not fundamental and do not need to be constant. Any slow change in the scale 
of the extra dimensions would be revealed by measurable changes in our three- 
dimensional 'constants'. Second, we appreciate that some apparent constant 
might be determined partially or completely by some spontaneous symmetry- 
breaking processes in the very early universe. This introduces an irreducible 
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random clement into the values of those constants. They may be different in 
different parts of the universe. The most dramatic manifestation of this process 
is provided by the chaotic and eternal inflationary universe scenarios. Third, 
any outcome of a theory of quantum gravity will be intrinsically probabilistic. 
It is often imagined that the probability distributions for observables will be 
very sharply peaked but this may not be the case for all possibilities. Thus, the 
value of G or G might be predicted to be spatially varying random variables. 
Fourth, the non-uniqueness of the vacuum state for the universe would allow 
other deals of the constants to have occurred in different places. At present we 
have no idea why any of the constants of Nature take the numerical values they 
do. Fifth, the observational limits on possible variations are often very weak 
(although they can be made to sound strong by judicious parametrisations). For 
example, the cosmological limits on varying G tell us only that G/G < 10~ 2 Hq, 
where Ho is the present Hubble rate. However, the last reason to consider 
varying constants is currently the most compelling. For the first time there is 
a body of detailed astronomical evidence for the time variation of a traditional 
constant. The observational programme of Webb et al @, [| has completed 
detailed analyses of three separate quasar absorption line data sets taken at 
Keck and finds persistent evidence consistent with the fine structure constant, 
a, having been smaller in the past, at z = 1 — 3.5. The shift in the value of a 
for all the data sets is given provisionally by Aa/a = (— 0.66± 0.11) x 10~ 5 . 
This result is currently the subject of detailed analysis and reanalysis by the 
observers in order to search for possible systematic biases in the astrophysical 
environment or in the laboratory determinations of the spectral lines. 

The first investigations of time-varying constants were those made by Lord 
Kelvin and others interested in possible time-variation of the speed of light at 
the end of the nineteenth century. In 1935 Milne devised a theory of gravity, 
of a form that we would now term 'bimetric', in which there were two times 
- one (t) for atomic phenomena, one (r) for gravitational phenomena - linked 
by t — log(t/io)- Milne Q required that the 'mass of the universe' (what we 
would now call the mass inside the particle horizon M w c?G~ l t) be constant. 
This required G oc t. Interestingly, in 1937 the biologist J.B.S. Haldane took a 
strong interest in this theory and wrote several papers || exploring its conse- 
quences for the evolution of life. The argued that biochemical activation energies 
might appear constant on the t timescale yet increase on the r timescale, giving 
rise to a non- uniformity in the evolutionary process. Also at this time there 
was widespread familiarity with the mysterious 'large numbers' O(10 40 ) and 
O(10 80 ) through the work of Eddington (although they had first been noticed 
by Weyl fl - see ref. and 0] for the history). These two ingredients were 
merged by Dirac in 1937 in a famous development (supposedly written on his 
honeymoon) that proposed that these large numbers (10 ) were actually equal, 
up to small dimcnsionlcss factors. Thus, if we form N ~ c 3 t/Gm n ~ 10 80 , 
the number of nucleons in the visible universe, and equate it to the square of 
N\ ~ e 2 /GmJ ~ 10 40 , the ratio of the electrostatic and gravitational forces 
between two protons then we are led to conclude that one of the constants, 
e,G,c,h,m n must vary with time. Dirac H chose G cc t^ 1 to carry the time 
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variation. Unfortunately, this hypothesis did not survive very long. Edward 
Teller || pointed out that such a steep increase in G to the past led to huge in- 
creases in the Earth's surface temperature in the past. The luminosity of the sun 
varies as L oc G 7 and the radius of the Earth's orbit as R oc G -1 so the Earth's 
surface temperature T e varies as (L/R 2 ) 1 / 4 oc G 9 / 4 oc t~ 9 / 4 and would exceed 
the boiling point of water in the pre-Cambrian era. Life would be eliminated. 
Gamow subsequently suggested that the time variation needed to reconcile the 
large number coincidences be carried by e rather than G, but again this strong 
variation was soon shown to be in conflict with geophysical and radioactive de- 
cay data. This chapter was brought to an end by Dickc JTo| who pointed out 
that the N ~ iVf large number coincidence was just the statement that t, the 
present age of the universe when our observations are being made, is of order 
the main sequence stellar lifetime, t ms ~ (Gm 2 l /hc)^ 1 h/m n c 2 ~ 10 10 yrs, and 
therefore inevitable for observers made from elements heavier than hydrogen 
and helium. Dirac never accepted this anthropic explanation for the large num- 
ber coincidences but curiously can be found making exactly the same type of 
anthropic argument to defend his own varying G theory by highly improbable 
arguments (that the Sun accretes material periodically during its orbit of the 
galaxy and this extra material cancels out the effects of overheating in the past) 
in correspondence with Gamow in 1967 (see JO for fuller details). 

Dirac's proposal acted as a stimulus to theorists, like Jordan, Brans and 
Dicke to develop rigorous theories which included the time variation of 
G self-consistently by modelling it as arising from the space-time variation of 
some scalar field c/>(x, t) whose motion both conserved energy and momentum 
and created its own gravitational field variations. In this respect the geometric 
structure of Einstein's equations provides a highly constrained environment to 
introduce variations of 'constants'. Whereas in Newtonian gravity we are at 
liberty to introduce a time-varying G(t) into the law of gravity by 

r z 

This creates a non-conservative dynamical system but can be solved fairly 
straightforwardly [Q. However, this strategy of simply 'writing in' the variation 
of G by merely replacing G by G(t) in the equations that hold when G is a 
constant fails in general relativity. If we were to imagine the Einstein equations 
would generalise to (G a b is the Einstein tensor) 

_ 8nG(t) 

W& — ~j + ab {*) 

then taking a covariant divergence and using V a G a t = 0, together with energy- 
momentum conservation (V a T a ;, = 0) requires that VG = and no variations 
are possible in eq. (^). Brans-Dicke theory is a familiar example of how the 
addition of an extra piece to T a b together with the dynamics of a G(cA) fields 
makes a varying G theory possible. Despite the simplicity of this lesson in 
the context of a varying G theory the lesson was not taken on board when 
considering the variations of other non-gravitational constants and the literature 
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is full of limits on their possible variation which have been derived by considering 
a theory in which the time- variation is just written into the equations which hold 
when the constant does not vary. Recently, the interest in the possibility that a 
varies in time has led to the first extensive exploration of simple self-consistent 
theories in which a variations occur through the variation of some scalar field. 



2 Brans-Dicke Theories 
2.1 Equations and solutions 



Consider the paradigmatic case of Brans-Dicke (BD) theory [jiTJ to fix theoret- 
ical ideas about varying G. The form of the general solutions to the Friedmann 
metric in BD theories are fully understood [TjJ, JlJ]. There are three essential 
field equations for the evolution of BD scalar field <j)(t) and the expansion scale 
factor a(t) in a BD universe 



87175 <f> i lubd <t> k 
a 4> 



* + 3-£ = —^—(P-3P) (4) 

p + i- a {p + P ) = (5) 

Here, u>bd is the BD constant parameter and the theory reduces to general 
relativity in the limit lubd — * 00 and (f> = G^ 1 —>■ constant. A general feature 
of the BD field equations is that any solution of general relativity for which the 
energy momentum tensor of matter has vanishing trace (eg vacuum, black body 
radiation, Yang-Mills, or magnetic field) is a particular (<fi = constant) solution 
of BD theory. 

The general solutions begin at high density dominated by the BD scalar field 
4> ~ G" 1 and approximated are well approximated by the spatially flat vacuum 
(p = p = 0) solutions: 



a(t) = ^/(A+D (6) 

m = <^ V( ^ (7) 

A = 1 + y/1 + 2u BD /Z (8) 

^BD 

This vacuum solution is the t — > attractor for the perfect-fluid solutions. The 
general perfect-fluid solutions with equation of state p = Tp and k = can all 
be found. At early times they approach the vacuum solutions but at late time 
they approach particular power-law exact solutions ]l5| : 

a(t) = t[ 2 + 2 "B-D( 1 - r )]/[ 4 +3"-BD(i-r 2 )] ^ 
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0(t) = ^ of P(i-3r)/[4 + 3 W (i-r=)] (10) 

At late times they approach particular exact power-law solutions for a(t) and 
(j)(t) and the evolution is 'Machian' in the sense that the cosmological evolution 
is driven by the matter content rather than by the kinetic energy of the free </> 
field. 

In the radiation era this particular solution is the standard general relativity 
solution: 

a(t) = t 1/2 ; (j>~ 1 oc G = constant (11) 

For p = the solutions have the form 

a(t) = t< 2 -™)/ 3 ; f'aGar, (12) 

which continues until the curvature term takes over the expansion. Here, n is 
related to the constant Brans-Dicke lobd parameter by 

n=z — I (13) 

4 + 3lu bd 

and the usual general relativistic Einstein de Sitter universe is obtained as 
u>bd — ► oo and n — > 0. 

In a curvature-dominated era of expansion, as t — > oo the solutions for 
T > —1/3 approach the general relativity Milne vacuum solution with 



a(t) = t (14) 
<j) oc G _1 = constant (15) 

Notice how the curvature domination turns off the variation of G. All pre- 
vious studies of varying G in cosmology have focussed on the k = models and 
have not noticed the important role that would be played by negative curvature. 
The existence of negative curvature can produce little residue of G variation in 
the universe today. It also highlights the usefulness of having constraints from 
different cosmic times and redshifts. 

The other important lesson to learn from the cosmological limits on varying 
G is that care must be taken when using local limits on lobd, say from light- 
bending by the Sun or the other solar system tests of general relativity, and 
then assuming that they can be used in cosmological models. In reality the 
evolution of the universe is inhomogeneous and there are very large variations 
in density between the solar system and the extragalactic universe. If we had a 
perfect numerical simulation of cosmology with a varying G wc would be able 
to determine the contours of G and G with position in the universe. Until we 
have more information of that sort from models it is unwise to assume that the 
rate of change of G in the solar system will be the same as it is on cosmological 
scales. 
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An interesting particular example of this problem is given by the power-law 
solutions above for the case with T = — 1. This is equivalent to the universe being 
dominated by a vacuum energy and leads to power-law accelerated expansion 
in BD theory with 

o(t) = t^ +u (16) 
G cc (p- 1 oc r 2 (17) 

Thus it appears that if our universe were to be expanding today (as observations 
of the recession of type I supernovae indicate (l6)) then G must be falling very 
rapidly (even faster than Dirac predicted) locally. Clearly this is impossible 
observationally jL7|. The flaw in the argument is that the p = —p stress does 
not dominate the dynamics of the solar system and we must not apply the 
cosmological solution for the variation of G in the solar system anymore than 
we should apply the cosmological solution for the variation of p to the solar 
system. 

It would be very interesting to find realistic solutions (exact, approximate, 
or numerical) for inhomogeneous cosmological models with 0(x, i) in order to 
obtain some perspective on the likely variation in the change in G from solar 
system to galactic and extragalactic scales. At present potentially the strongest 
cosmological limit on time- varying G is stronger than solar system tests and 
comes, somewhat surprisingly, from the power spectrum of galaxy clustering. 
The effect of varying G is to shift the cosmic epoch of equality between the 
matter and radiation densities which determines the location of the peak of the 
clustering power spectrum JL8[ . 

3 A Simple Varying-Alpha Theory 



We are going to consider some of the cosmological consequences of a simple 
theory with time varying a. Such a theory was first formulated by Bekenstein 
[ fl9| as a generalisation of Maxwell's equations but ignoring the consequences 
for the gravitational field equations. Recently, Magueijo, Sandvik and myself 
have completed this theory |2^, ^l], Q to include the coupling to the 
gravitational sector and analysed its general cosmological consequences. This 
theory considers only a variation of the electromagnetic coupling and so far 
ignores any unification with the strong and electroweak interactions. We shall 
not discuss simultaneous variation of the electromagnetic and gravitational con- 
stants although that analysis can be done and is presented elsewhere (during 
the dust era of a flat Friedmann universe with varying a(t) and G(i), their time- 
evolution approaches an attractor in which the product aG is a constant and 
a on G^ 1 oc t n , where n is given by eq. (|l3|).) 

Our aim in studying this theory is to build up understanding of the effects of 
the expansion on varying a and to identify features that might carry over into 
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more general theories in which all the unified interactions vary |25], g(| g7[ . The 
constraint imposed on varying a by the need to bring about unification at high 
energy is likely to be significant but the complexities of analysing the simulta- 
neous variation of all the constants involved in the supersymmetric version of 
the standard model are considerable. At the most basic level we recognise that 
any time variation in the fine structure could be carried by either or both of the 
electromagnetic or weak couplings above the electroweak scale. 

The idea that the charge on the electron, or the fine structure constant, might 
vary in cosmological time was proposed in 1948 by Teller, ||, who suggested that 
a oc (lni) -1 was implied by Dirac's proposal that G oc t" 1 and the numerical 
coincidence that a -1 ~ hx{hc/ Gm^ r ) , where the proton mass. Later, in 

1967, Gamow p^] suggested a oc t as an alternative to Dirac's time-variation 
of the gravitation constant, G, as a solution of the large numbers coincidences 
problem and in 1963 Stanyukovich had also considered varying a, p9| , in this 
context. However, this power-law variation in the recent geological past was 
soon ruled out by other evidence |?0| . 

There are a number of possible theories allowing for the variation of the 
fine structure constant, a. In the simplest cases one takes c and H to be con- 
stants and attributes variations in a to changes in e or the permittivity of free 
space (see Q for a discussion of the meaning of this choice). This is done 
by letting e take on the value of a real scalar field which varies in space and 
time (for more complicated cases, resorting to complex fields undergoing spon- 
taneous symmetry breaking, see the case of fast tracks discussed in p2[). Thus 
eo — > e = eoe(a; M ), where e is a dimensionless scalar field and eo is a constant 
denoting the present value of e. This operation implies that some well estab- 
lished assumptions, like charge conservation, must give way ]35| ]. Nevertheless, 
the principles of local gauge invariance and causality are maintained, as is the 
scale invariance of the e field (under a suitable choice of dynamics). In addition 
there is no conflict with local Lorentz invariance or covariance. 

With this set up in mind, the dynamics of our theory is then constructed as 
follows. Since e is the electromagnetic coupling, the e field couples to the gauge 
field as eA^ in the Lagrangian and the gauge transformation which leaves the 
action invariant is eA^ — ► eA^ +X ui rather than the usual A^ — ► A^ +X,a- The 
gauge-invariant electromagnetic field tensor is therefore 

= \ ((eA v ) >fl - (eAJ, v ) , (18) 

which reduces to the usual form when e is constant. The electromagnetic part 
of the action is still 

S em = -J d 4 xV=gF^F^. (19) 

and the dynamics of the e field are controlled by the kinetic term 

(20) 
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as in dilaton theories. Here, I is the characteristic length scale of the theory, 
introduced for dimensional reasons. This constant length scale gives the scale 
down to which the electric field around a point charge is accurately Coulombic. 
The corresponding energy scale, Hc/l, has to lie between a few tens of MeV and 
Planck scale, ~ l0 19 GeV to avoid conflict with experiment. 

Our generalisation of the scalar theory proposed by Bekenstein [|l9j described 
in ref. |^2|, Sll ^0], includes the gravitational effects of ij) and gives the field 
equations: 

Gfiv = SirG (T™ a " er + It + T ™^) ■ (21) 

The stress tensor of the ip field is derived from the lagrangian C$ = —^d^ipd^ip 
and the ip field obeys the equation of motion 

Op = -e-^£ em (22) 

where we have defined the coupling constant to = {c)/l 2 . This constant is of 
order ~ 1 if, as in [?], the energy scale is similar to Planck scale. It is clear 
that C em vanishes for a sea of pure radiation since then C em — (E 2 — B 2 )/2 = 
0. We therefore expect the variation in a to be driven by electrostatic and 
magnetostatic energy-components rather than electromagnetic radiation. 

In order to make quantitative predictions we need to know how much of the 
non-relativistic matter contributes to the RHS of Eqn. (p2|). This is parametrised 
by £ = Cem I p, where p is the energy density, and for baryonic matter C em = 
E 2 /2. For protons and neutrons C p an d Cn can be estimated from the electro- 
magnetic corrections to the nucleon mass, 0.63 MeV and —0.13 MeV, respec- 
tively Jt2| . This correction contains the E 2 /2 contribution (always positive), 
but also terms of the form j^a^ (where j M is the quarks' current) and so cannot 
be used directly. Hence we take a guiding value Q p ( n ~ 10 -4 . Furthermore 
the cosmological value of £ (denoted £„J has to be weighted by the fraction of 
matter that is non-baryonic, a point ignored in the literature ]l9| . Hence, £ m 
depends strongly on the nature of the dark matter and can take both positive 
and negative values depending on which of Coulomb-energy or magnetostatic 
energy dominates the dark matter of the Universe. It could be that Ccdm ~ — 1 
(superconducting cosmic strings, for which C em s» —B 2 /2), or Ccdm ^ 1 (neu- 
trinos). BBN predicts an approximate value for the baryon density of ss 0.03 
with a Hubble parameter of ho ~ 0.6 , implying £Icdm ~ 0.3. Thus depending 
on the nature of the dark matter C m can be virtually anything between — 1 and 
+ 1. The uncertainties in the underlying quark physics and especially the con- 
stituents of the dark matter make it difficult to impose more certain bounds on 

Cm " 

We should not confuse this theory with other similar variations. Bekenstein's 
theory does not take into account the stress energy tensor of the dielectric field 
in Einstein's equations, and their application to cosmology. Dilaton theories 
predict a global coupling between the scalar and all other matter fields. As a 
result they predict variations in other constants of nature, and also a different 
dynamics to all the matter coupled to electromagnetism. An interesting appli- 
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cation of our approach has also recently been made to braneworld cosmology in 



3.1 The cosmological equations 

Assuming a homogeneous and isotropic Friedmann metric with expansion scale 
factor a(t) and curvature parameter k in eqn. (pl|), we obtain the field equations 
(c=l) 



2 

a ' 



3 

k A 
^ + 3 



(pm (1 + Cm exp [-2^]) + P r exp [-2$ + 



(23) 



where A is the cosmological constant. For the scalar field we have the propaga- 
tion equation, 

$ + 3Hi>=--exp[-2il>K m P m , (24) 

where H = a/a is the Hubble expansion rate. We can rewrite this more simply 
as 

(>a 3 ) = Nexp[-2V>] (25) 
where N is a positive constant defined by 

jV = _ 2 CmPm° 3 (26) 

Note that the sign of the evolution of -0 is dependent on the sign of £ m . 
Since the observational data is consistent with a smaller value of a in the past, 
we will in this paper confine our study to negative values of £ m , in line with our 
recent discussion in Refs. ]2^, |2(], The conservation equations for the 
non-interacting radiation and matter densities are 

p m + 3Hp m = (27) 
p r +4Hp r = 2i>p r . (28) 

and so p m cx a~ 3 and p r e~ 2 ^ cx a -4 , respectively. If additional non-interacting 
perfect fluids satisfying equation of state p = (j—l)p are added to the universe 
then they contribute density terms p cx a -37 to the RHS of eq.fl23|) as usual. 
This theory enables the cosmological consequences of varying e, to be analysed 
self-consistently rather than by changing the constant value of e in the standard 
theory to another constant value, as in the original proposals made in response 
to the large numbers coincidences. 

We have been unable to solve these equations in general except for a few 
special cases J35| . However, as with the Friedmann equation of general relativity, 
it is possible to determine the overall pattern of cosmological evolution in the 
presence of matter, radiation, curvature, and positive cosmological constant by 
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matched approximations. We shall consider the form of the solutions to these 
equations when the universe is successively dominated by the kinetic energy of 
the scalar field ip, pressure-free matter, radiation, negative spatial curvature, 
and positive cosmological constant. Our analytic expressions are checked by 
numerical solutions of (S3) and (pi[). 

3.2 Observational implications 

There are a number of conclusions that can be drawn from the study of the 
simple BSBM models with £ TO < 0. These models give a good fit to the varying 
a implied by the QSO data of refs. B|. There is just a single parameter to 
fit and this is given by the choice 

= (2±1) x 1(T 4 (29) 

UJ 

The simple solutions predict a slow (logarithmic) time increase during the 
dust era of k — Friedmann universes. The cosmological constant turns off 
the time-variation of a at the redshift when the universe begins to accelerate 
(z ~ 0.7) and so there is no conflict between the a variation seen in quasars at 
z ~ 1 — 3.5 and the limits on possible variation of a deduced from the operation 
of the Oklo natural reactor |3(| (even assuming that the cosmological variation 
applies unchanged to the terrestrial environment). The reactor operated 1.8 
billion years ago at a redshift of only z ~ 0.1 when no significant variations 
were occurring in a. The slow logarithmic increase in a also means that we 
would not expect to have seen any effect yet in the anisotropy of the microwave 
backgrounds |57| pSf : the value of a at the last scattering redshift, z — 1000, 
is only 0.005% lower than its value today, similarly, the essentially constant 
evolution of a predicted during the radiation era leads us to expect no measur- 
able effects on the products of Big Bang nucleosynthesis (BBN) J3?| because a 
was only 0.007% smaller at BBN than it is today. This does not rule out the 
possibility that unification effects in a more general theory might require varia- 
tions in weak and strong couplings, or their contributions to the neutron-proton 
mass difference, which might produce observable differences in the light element 
productions and new constraints on varying a at z ~ 10 9 — 10 10 . By contrast 
varying alpha cosmologies with £ > lead to bad consequences . The fine struc- 
ture falls rapidly at late times and the variation is such that it even comes to 
dominate the Friedmann equation for the cosmological dynamics. We regard 
this as a signal that such models are astrophysically ruled out and perhaps also 
mathematically badly behaved. 

We should also mention that theories in which a varies will in general lead 
to violations of the weak equivalence principle (WEP). This is because the a 
variation is carried by a field like ip and this couples differently to different 
nuclei because they contain different numbers of electrically charged particles 
(protons). The theory discussed here has the interesting consequence of leading 
to a relative acceleration of order 10" 13 [fio) if the free coupling parameter is 
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fixed to the value given in eq. (^9|) using a best fit of the theories cosmological 
model to the QSO observations of refs. Jjl. Other predictions of such WEP 
violations have also been made in refs. [[Il], [h| Q. The observational 
upper bound on this parameter is just an order of magnitude larger, at 10 -12 , 
but space-based tests planned for the STEP mission are expected to achieve a 
sensitivity of order 10 -18 and will provide a completely independent check on 
theories of time- varying e and a. This is an exciting prospect for the future. 

3.3 The nature of the Priedmann solutions 

The cosmological behaviour of the solutions to these equations was studied 
by us in detail, both analytically and numerically in refs. |22|, EH], |2^, p4[ , 
p5) . Typically, the variation in a does not have a significant effect on the 
evolution of the scale factor at late times although the cosmological expansion 
does significantly affect the evolution of a. The evolution of a is summarised as 
follows: 

During the radiation era a(t) ~ t 1 / 2 and a is constant in universes with our 
entropy per baryon and present value of a like our own. It increases in the dust 
era, where a(t) ~ i 2 / 3 . The increase in a however, is very slow with a late-time 
solution for if> proportional to \ log(2iVlog(i)), and so 

a ~2N log t (30) 

This slow increase continues until the expansion becomes dominated by neg- 
ative curvature, a(t) ~ t, or by a cosmological vacuum energy, a(t) ~ exp[Af/3]. 
Thereafter a asymptotes rapidly to a constant. If we set the cosmological con- 
stant equal to zero and k — then, during the dust era, a would continue to 
increase indefinitely. The effect of the expansion is very significant at all times. 
If we were to turn it off and set a(t) constant then we could solve the ip equation 
to give the following exponentially growing evolution for a, f35j| : 

a = exp[2?/>] = A~ 2 cosh 2 L4iV 1/2 (i + t )}; A constant. (31) 

From these results it is evident that non-zero curvature or cosmological con- 
stant brings to an end the increase in the value of a that occurs during the 
dust-dominated era. Hence, if the spatial curvature and A are both too small 
it is possible for the fine structure constant to grow too large for biologically 
important atoms and nuclei to exist in the universe. There will be a time in the 
future when a reaches too large a value for life to emerge or persist. The closer 
a universe is to flatness or the closer A is to zero so the longer the monotonic 
increase in a will continue, and the more likely it becomes that life will be ex- 
tinguished. Conversely, a non-zero positive A or a non-zero negative curvature 
will stop the increase of a earlier and allow life to persist for longer. If life can 
survive into the curvature or A-dominated phases of the universe's history then 
it will not be threatened by the steady cosmological increase in a unless the 
universe collapses back to high density. 
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This type of behaviour can also be found in the presence of time-varying 
G. If a BD dust universe is exactly flat (fc = 0) then G will continue to fall 
forever. Only if there is negative curvature will the evolution of G eventually be 
turned off and the expansion asymptote to the Milne behaviour with a = t and 
G — > constant. Again, without the small deviation from flatness the strength of 
gravity would ultimately become too weak for the existence of stars and planets 
and the universe would become biologically inhospitable, if not uninhabitable. 

There have been several studies, following Carter, |Q and Tryon of 
the need for life-supporting universes to expand close to the 'flat' Einstein de 
Sitter trajectory for long periods of time. This ensures that the universe cannot 
collapse back to high density before galaxies, stars, and biochemical elements can 
form by gravitational instability, or expand too fast for stars and galaxies to form 
by gravitational instability fiq , ffl]. Likewise, it was pointed out by Barrow and 
Tipler, Q that there are similar anthropic restrictions on the magnitude of any 
cosmological constant, A. If it is too large in magnitude it will cither precipitate 
premature collapse back to high density (if A < 0) or prevent the gravitational 
condensation of any stars and galaxies (if A > 0). Thus existing studies provide 
anthropic reasons why we can expect to live in an old universe that is neither 
too far from flatness nor dominated by a much stronger cosmological constant 
than observed (|A| < 10 |A obs |). 

Inflationary universe models provide a possible theoretical explanation for 
proximity to flatness but no explanation for the smallncss of the cosmologi- 
cal constant. Varying speed of light theories (4^, |3l], |o[ ||l], [52| offer possible 
explanations for proximity to flatness and smallness of a classical cosmologi- 
cal constant (but not necessarily for one induced by vacuum corrections in the 
early universe). We have shown that if we enlarge our cosmological theory to 
accommodate variations in some traditional constants then it appears to be an- 
thropically disadvantageous for a universe to lie too close to flatness or for the 
cosmological constant to be too close to zero. This conclusion arises because of 
the coupling between time-variations in constants like a and the curvature or 
A, which control the expansion of the universe. The onset of a period of A or 
curvature domination has the property of dynamically stabilising the constants, 
thereby creating favourable conditions for the emergence of structures. This 
point has been missed in previous studies because they have never combined 
the issues of A and flatness and the issue of the values of constants. By coupling 
these two types of anthropic considerations we find that too small a value of A 
or the spatial curvature can be as poisonous for life as too much. Universes like 
those described above, with increasing a(t), lead inexorably to an epoch where 
a is too large for the existence of atoms, molecules, and stars to be possible. 

Surprisingly, there has been almost no consideration of habitability in cos- 
mologies with time- varying constants since Haldane's discussions Q of the bio- 
logical consequences of Milne's bimetric theory of gravity. Since then, attention 
has focussed upon the consequences of universes in which the constants are dif- 
ferent but still constant. Those cosmologies with varying constants that have 
been studied have not considered the effects of curvature or A domination on 
the variation of constants and have generally considered power-law variation to 
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hold for all times. The examples described here show that this restriction has 
prevented a full appreciation of the coupling between the expansion dynamics 
of the universe and the values of the constants that define the course of local 
physical processes within it. Our discussion of a theory with varying a shows 
for the first time a possible reason why the 3-curvature of universes and the 
value of any cosmological constant may need to be bounded below in order that 
the universe permit atomic life to exist for a significant period. Previous an- 
thropic arguments have shown that the spatial curvature of the universe and 
the value of the cosmological constant must be bounded above in order for life- 
supporting environments (stars) to develop. We note that the lower bounds 
discussed here are more fundamental than these upper bounds because they de- 
rive from changes in a which have direct consequences for biochemistry whereas 
the upper bounds just constrain the formation of astrophysical environments by 
gravitational instability. Taken together, these arguments suggest that within 
an ensemble of all possible worlds where a and G are time variables, there might 
only be a finite interval of non-zero values of the curvature and cosmological 
constant contributions to the dynamics that both allow galaxies and stars to 
form and their biochemical products to persist. 

3.4 The role of inhomogeneities 

We can also detect where and how we might expect spatial variations to arise 
in a fuller description. Aside from the complexities of the full inhomogeneous 
cosmological solution for the formation of galaxies, stars, and planets, we can 
isolate non-uniformities that enter through the constant parameter N which 
dictates the form and time-evolution of a(t) and a(t). First we see that N is 
proportional to the density of electromagnetically charged matter in the uni- 
verse. This will possess some spatial variation and is of order on large 
scales. More significant though is the variation of the baryonic content of the 
CDM density with scale. We need the CDM to be dominated by matter with 
magnetic charge (but see Bekenstein). This can be the case on large scale but 
we know that the dark matter becomes dominated by baryons (therefore with 
C > 0) locally. Hence, there is expected to be a very significant spatial variation 
of £ with scale, including a change of sign, which will feed into the variation of 
a. 

3.5 General properties of the evolution of alpha and G 

The evolution equation for tp(t) has a number of simple but important proper- 
ties. Since N > the right-hand side or eq. ( p5| ) must be positive. This means 
that linearisations of this equation are dangerous and give rise to linearisation 
instabilities unless attention is confined to the regime ip « 1. In general the 
positivity property means that there can be no oscillations of ip or a in time in 
solutions of this equation. This follows from the required positivity of (tpa 3 ), 
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which means that ^cannot have a maximum. The possible cosmological evolu- 
tions for ip and a are decrease to a minimum followed by a monotonic increase, 
monotonic decrease, or monotonic increase. This conclusion holds independently 
of the value of k in the Friedmann equation. This has one very important con- 
sequence. It means that the asymptotic monotonic non-decrease of a found in 
our flat and open universes will still occur in closed universes. There cannot 
be a sudden change in the evolution of a when the universe starts to collapse. 
This also means that if we model spherical overdensities by closed universes 
embedded in a flat background then the evolution of a(t) in the overdensities 
will be very similar to that in the flat background even when the overdensities 
collapse to form bound 'clusters'. This has the important implication that such 
an inhomogeneous universe will not end up with very different values of a and 
a in inside and outside the bound inhomogeneities. 

This argument can also be applied to the evolution of <f> and G in BD theory. 
Consider the case of dust (p = 0). The combination (<f)a 3 ) must now be positive 
and so <f> cannot have a maximum and G cannot have a minimum regardless 
of the sign of the curvature parameter k. In particular, G(t) cannot oscillate. 
Again, this property acts as a safeguard on the divergent evolution of G inside 
and outside overdensities. 



4 The Second Law 

There has been considerable recent discussion |53|, Q about the equivalence 
of models of the variation of different dimensional 'constants' of Nature. In 
particular, it has been suggested that consideration of the second law of black 
hole thermodynamics distinguishes, say, variations of e from variations of c and 
that some of these variations could be ruled out because they bring about a 
decrease in time of the Bckenstein-Hawking entropy of a charged black hole. 
Others have argued that no such distinction is operationally possible. However, 
we believe that the most crucial factor has been missed in this discussion. In 
theories which generalise general relativity by allowing traditional constants 
(like G or e) to vary the black hole solutions with event horizons are particular 
solutions of the theory in which the constant concerned is a constant. When the 
constant varies the black hole solution no longer exists and there is no longer 
any black hole thermodynamics to constrain the variation. The situation is very 
clear in the simple case of a Schwarzschild black hole in Brans Dicke theory. We 
know from the work of Hawking |55| that the black hole solutions are the same as 
those in general relativity. Thus Schwarzschild is a 0~G _1 = constant solution 
of the Brans-Dicke field equations. The entropy of this black hole is 

Sbh GM 2 . 

If we were to apply the second law to this formula it would appear to say that all 
cosmological solutions in which G falls with time are ruled out. However, this 
would not be a correct deduction (which is fortunate because we see from eq. 
(fl0|) that essentially all Brans-Dicke cosmologies have such behaviour) because 
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<f> and G are constant on the Schwarzschild horizon. If we allow variation of 
G then the solution turns into a naked singularity and the thermodynamic 
relations no longer exist. Thus one cannot at present use considerations of black 
hole thermodynamics to constrain or distinguish the time or space variation of 
constants of Nature by simply 'writing in' time variations into the formulae that 
define the black hole when these constants do not vary. 
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